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CLAIMS 

1 . (Original) A method of digitally processing data in a data array defining a target matrix (X) 
using non-negative matrix factorisation to determine a pair of matrices (F, G), a first matrix of 
said pair determining a set of features for representing said data, a second matrix of said pair 
determining weights of said features, such that a product of said first and second matrices 
approximates said target matrix, the method comprising: 

inputting said target matrix data (X); 

selecting a row of said one of said first and second matrices and a column of the other of 
said first and second matrices; 

determining a target contribution (R) of said selected row and column to said target 

matrix; 

determining, subject to a non-negativity constraint, updated values for said selected row 
and column from said target contribution; and 

repeating said selecting and determining for the other rows and columns of said first and 
second matrices until all said rows and columns have been updated. 

2. (Original) A method as claimed in claim 1 wherein said determining of updated values 
comprises determining a new value for said selected row substantially independent of a previous 
value of said selected row and determining a new value for said selected column substantially 
independent of a previous value of said selected column. 

3. (Previously Presented) A method as claimed in claim 1, wherein said determining of a target 
contribution (R) comprises determining a difference between said target matrix (X) and a sum of 
weighted features determined from all said rows and columns of said first and second matrices 
except said selected row and column. 

4. (Previously Presented) A method as claimed in claim 1 , wherein said determining of updated 
values comprises calculating values of Gj a = fi (R,F) and F au = h (R, G) where R is a matrix with 
I rows and U columns, F is a matrix with A rows and U columns and G is a matrix with I rows 
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and A columns, where fi and f~2 denote first and second functions, and where Gj a denotes a data 
element in the zth row and ath column of G and F au denotes a data element in the ath row and wth 
column of F. 



5. (Original) A method as claimed in claim 4 wherein fi and f2 are selected to minimise a cost 
function measuring a quality of approximation of a product of said selected row and column to 
said target contribution. 

6. (Original) A method as claimed in claim 5 wherein said cost function comprises a squared 
Euclidean distance between a product of said selected row and column and said target 
contribution. 



7. (Original) A method as claimed in claim 5 wherein said cost function comprises a divergence 
function between a product of said selected row and column and said target contribution. 

8. (Previously Presented) A method as claimed in claim 4, wherein Gj a and Fj U are determined in 
accordance with: 

U I 

Yr. f o. Yg. r. o. 

/ j iu au m / j la m iu 

(l = ujzl 77 -j^l 

w za U » au I 

Vf 2 0. Yg 2 0. 

/ j au iu / j m iu 

u = \ 1=1 

where R; u denotes a data element in the zth row and uth column of R, where R is given by: 

A,n*a 

R IU = X IU ~ X G In F nU 

n = l 

and where O iu denotes a data element in the zth row and z4h column of an I by U matrix O. 



9. (Original) A method as claimed in claim 8 wherein Oiu is substantially unity for all i and u. 
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10. (Original) A method as claimed in claim 8 wherein <3>i u has the form <Dj u = I/(Z; U + y) where 
Zi U denotes a data element in the z'th row and uth column of an I by U matrix dependent upon at 
least one of X and R, and y is positive. 

1 1 . (Previously Presented) A method as claimed in claim 4, wherein A is less than the smaller 
of I and U. 

12. (Previously Presented) A method as claimed in claim 1, further comprising initialising said 
first and second matrices. 

13. (Original) A method as claimed in claim 12 wherein said data comprises image data for an 
image in a time series of images, and wherein said initialising is conditional upon a degree of 
difference between said image and a previous image. 

14. (Previously Presented) A method as claimed in claim 1, wherein said determining subject to 
a non-negativity constraint comprises setting a said updated value to substantially zero where the 
updated value would otherwise be negative. 

15. (Previously Presented) A method as claimed in claim 1, further comprising constraining said 
updated values to lie between a minimum and a maximum value. 

16. (Previously Presented) A method as claimed in claim 1, further comprising repeating said 
updating of all said rows and columns of said first and second matrices for a plurality of 
iterations. 

17. (Previously Presented) A method as claimed in claim 1, wherein said data comprises image 
data defining an image, and wherein said set of features determined by said first matrix 
comprises a set of subframes which when combined according to said weights determined by 
said second matrix approximate said image. 
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18. (Previously Presented) A method of driving a display comprising a plurality of pixels 
arranged in rows and columns, the method comprising employing the method of claim 1 to 
process data for display as said target matrix data (X) to determine said first and second matrices 
(F,G), and driving said display to form an image using a plurality of subframes, each subframe 
having said rows and columns of pixels driven responsive to a row of one of said first and second 
matrices and a column of the other of said first and second matrices. 

19. (Previously Presented) A method of image matching, the method comprising inputting 
target image data for matching as said target matrix data; processing said image data as claimed 
in claim 1 to determine said first and second matrices; comparing data from at least one of said 
first and second matrices with stored data comprising data for a corresponding first and/or second 
matrix for a second image; and outputting image match data responsive to said comparing. 

20. (Previously Presented) A method as claimed in claim 19 wherein said stored data comprises 
data for a plurality of images held in a non- volatile store in the form of at least one of said first 
matrix said second matrix. 

21. (Original) A method as claimed in claim 19 wherein said stored data comprises stored data 
derived from said second image, the method further comprising inputting said second image data 
and processing said data input to determine said first and second matrices for said second image. 

22. (Previously Presented) A method as claimed in claim 19, wherein both said target image 
data and said second image data define a biometric image, in particular a facial image. 

23. (Previously Presented) A method of providing security data for a data network, the method 
comprising: 

reading traffic on said network to provide data for analysis; 

processing said data for analysis as claimed in claim 1 to determine said first and second 
matrices; and 
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analysing data from at least one of said first and second matrices to determine security 
data for the network. 

24. (Original) A method as claimed in claim 23 further comprising controlling access to and/or 
traffic on said network responsive to said analysing. 

25. (Previously Presented) A method of data mining, the method comprising: 

applying the method of claim 1 to data stored in a database to determine a set of 
discovered features; and 

outputting analysis data derived from expression of said stored data in terms of said 
discovered features. 

26. (Previously Presented) A method of processing sensor data, the method comprising: 

inputting said sensor data as a data array; and 

processing said sensor data using the method of claim 1 ; and 

outputting said processed data expressed in terms of a set of said features. 

27. (Previously Presented) A method of analysing biological data, the method comprising: 

inputting said biological data; 

processing said biological data using the method of claim 1 to determine feature data 
representing features for said biological data; and 

analysing said biological data using said feature data. 

28. (Original) A method as claimed in claim 27 wherein said biological data comprises atomic 
coordinate data and wherein said features correspond to physically similar or complementary 
features of biological entities. 

29. (Original) A method as claimed in claim 27 wherein said biological data comprises 
sequence data in particular one or more of genome sequence data, proteome sequence data, 
microarray data, amino acid sequence data, and nucleotide sequence data, the method further 
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comprising outputting analysing data for matching portions of said biological data. 

30. (Previously Presented) A method of teaching a data processing system, the method 
comprising: 

inputting data object data for a plurality of instances of data objects about which said data 
processing system is to learn; 

processing said data object data using the method of claim 1 to identify feature data 
defining one or more features defining characteristics of said objects; and 

updating an information store of said data processing system using said feature data. 

31. (Original) A method as claimed in claim 30 wherein a number or dimension of said features 
is less than a number of attributes or dimension of said objects. 

32. (Previously Presented) A method of data analysis, the method comprising: 

inputting data for analysis; 

processing said data for analysis using the method of claim 1 to determine feature data 
representing a plurality of features of said data for analysis; and 

analysing said data for analysis by analysing said feature data. 

33. (Previously Presented) A carrier medium carrying processor control code, to, when running, 
implement the method of claim 1 . 

34. (Previously Presented) A computer system inputting said target matrix data (X), the system 
comprising: 

an input for said data for said data array; 

an output for outputting said first and second matrices; 

data memory for storing said target matrix and said pair of matrices; 

program memory storing processor control code; and 
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a processor coupled to said input, to said output, to said data memory and to said program 
memory for loading and implementing said processor control code, said code comprising code to, 
when running implement the method of claim 1 . 

35. (Original) Apparatus for digitally processing data in a data array defining a target matrix (X) 
using non-negative matrix factorisation to determine a pair of matrices (F, G), a first matrix of 
said pair determining a set of features for representing said data, a second matrix of said pair 
determining weights of said features, such that a product of said first and second matrices 
approximates said target matrix, the apparatus comprising: 

means for inputting said target matrix data (X); 

means for selecting a row of said one of said first and second matrices and a column of 
the other of said first and second matrices; 

means for determining a target contribution (R) of said selected row and column to said 
target matrix; 

means for determining, subject to a non-negativity constraint, updated values for said 
selected row and column from said target contribution; and 

means for repeating said selecting and determining for the other rows and columns of said 
first and second matrices until all said rows and columns have been updated. 

36. (Original) A method of driving an electro-optic display, the display having a matrix of 
pixels, the method comprising: 

inputting image data for said matrix of pixels into an image data matrix; 

factorising said image data matrix into a product of first and second factor matrices; and 

driving said display using said factor matrices; and wherein 

said factorising comprising iteratively adjusting said factor matrices such that their 
product approaches said image data matrix; and wherein 

said iterative adjusting comprises adjusting each row of one of said factor matrices and 
each column of the other of said factor matrices in turn. 
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37. (Original) A method as claimed in claim 36 wherein said adjusting of a row or column 
comprises determining a new value for a said row or column substantially independently of a 
previous value of said row or column. 

38. (Previously Presented) A method as claimed in claim 36, wherein each said pixel is 
addressable by a row electrode and a column electrode; wherein one of said factor matrices 
defines row drive signals and the other of said factor matrices defines column drive signals; and 
wherein said driving comprises driving using said row and column drive signals. 

39. (Original) A method as claimed in claim 38 wherein said driving comprises driving a 
plurality of said row electrodes in combination with a plurality of said column electrodes. 

40. (Previously Presented) A method as claimed in claim 38, wherein said driving comprises 
driving said display with successive sets of said row and column signals to build up a display 
image, each said set of signals defining a subframe of said display image, said subframes 
combining to define said display image. 

41 . (Original) A method as claimed in claim 40 wherein said first factor matrix has dimensions 
determined by a number of said row electrodes and a number of said subframes, and wherein said 
second factor matrix has dimensions determined by a number of said column electrodes and said 
number of subframes. 

42. (Previously Presented) A method as claimed in claim 40, wherein a number of said 
subframes is no greater than the smaller of a number of said row electrodes and a number of said 
column electrodes. 

43. (Previously Presented) A method as claimed in claim 36, wherein said display comprises a 
passive matrix OLED display. 
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44. (Previously Presented) A method as claimed in claim 36, wherein said display comprises a 
plasma display. 

45. (Previously Presented) A method as claimed in claim 36, wherein said display comprises an 
inorganic LED display. 

46. (Previously Presented) A method as claimed in claim 36, wherein said display comprises a 
liquid crystal display. 

47. (Previously Presented) A carrier medium carrying processor control code, to, when running, 
implement the method of claim 36 . 

48. (Original) A driver for an electro-optic display, the display having a matrix of pixels, the 
driver comprising: 

an input to input image data for said matrix of pixels into an image data matrix; 

a matrix factorisation system to factorise said image data matrix into a product of first 
and second factor matrices; and 

a driver output to drive said display using said factor matrices, and wherein 

said matrix factorisation system is configured to iteratively adjust said factor matrices, 
such that their product approaches said image data matrix, by adjusting each row of one of said 
factor matrices and each column of the other of said factor matrices in turn. 

49. (Original) An integrated circuit including the matrix factorisation system of claim 48. 

50. (Original) A method of processing a data array defining a target matrix (X) to determine a 
pair of factor matrices (F,G) such that a product of said factor matrices approximates said target 
matrix (X), the method comprising: 

determining for a single row or column of a first said factor matrix a value to which an 
updating rule would converge when iteratively applied, said updating rule comprising an 
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updating rule of a factorising algorithm which iteratively updates two factor matrices to more 
closely approximate a target matrix using said updating rule; 

updating said row or column with said determined value; 

repeating said determining and updating for a column or row of a second said factor 
matrix; and 

repeating said determining and updating of said first and second factor matrices to update 
each row or column of said first factor matrix and each column or row of said second factor 
matrix. 

5 1 . (Original) A carrier carrying processor control code to, when running, implement the 
method of claim 50. 



